BCS theory for s+g-wave superconductivity in borocarbides Y(Lu)Ni2B2C 
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The s+g mixed gap function Ak = A[(l — x) — x sin* cos 4;/!)] (x: weight of g-wave component) 
has been studied within BCS theory. By suitable consideration of the pairing interaction, we have 
confirmed that the coexistence of s- and g-wave, as well as the state with equal s and g amplitudes 
(i.e., X — 1/2) may be stable. This provides the semi-phenomenological theory for the s+g-wave su- 
perconductivity with point nodes which has been observed experimentally in borocarbides YNi2B2C 
and possibly in LuNi2B2C. 



PACS numbers: 74.20.Rp, 74.20.Fg, 74.70. Dd 
I. INTRODUCTION 

The rare earth nickel borocarbides RNi2B2C (R=Y, 
Lu, Tm, Er, Ho, and Dy) have attracted great interest 
in recent years due to superconductivity (SC) as well as 
its possible coexistence with antiferromagnetic order. ^'^ 
It has initially been thought that these materials can be 
understood by an isotropic s-wave pairing via the conven- 
tional electron-phonon coupling."^ However, recent var- 
ious experimental results particularly on the two non- 
magnetic borocarbides Y(Lu)Ni2B2C, including specific 
heat,^~^ thermal conductivity,^'^ Raman scattering, 
NMR relaxation rate,^^ photoemission spectroscopy,^^ 
scanning tunneling microscopy and spectroscopy,^^ have 
unambiguously shown that the gap function is highly 
anisotropic with anisotropy ratio^ Amin/^max <10^^. 
For example, the ^/H dependence of the specific heat in 
the vortex state indicates a superconducting state with 
nodal excitations''"®. The power law behavior of the 
spin-lattice relaxation rate^^ also suggests the existence 
of nodes. Very recently, compelling evidence is presented 
by Izawa et al. from the angular-dependent thermal 
conductivity in a magnetic field that the gap function 
of YNi2B2C has point nodes which are located along 
the [1,0,0] and [0,1,0] directions.^ The same conclusion 
can be also drawn from the angular-dependent specific 
heat data.^ Highly anisotropic s-wave gap (with possi- 
ble nodes) was also discovered in LuNi2B2C by thermal 
conductivity measurements as a function of temperature 
and field strength.^ Thus the previous isotropic s-wave 
theory has to be critically reconsidered. 

Recently Maki et al. have proposed that the so called 
s-|-g-wave spin singlet gap function for Y(Lu)Ni2B2C su- 
perconductors, i.e.,^^'^^'^ 

Ak = y(l-sin4^?cos4,/)) , (1) 

is consistent with experimental results. We introduced 
9, cj) as polar and azimuthal angles of k, respectively. 
Here the second 'g-wave' contribution is given by a fourth 
degree fully symmetric (Aig) basis function ■ip^'^\9,4>) 
in tetragonal D4/1 symmetry which is, up to a con- 
stant, equal to the real 'tesseral harmonic' Z!l^{9,4>) = 
{l/V2)[Y^{9,(j)) +Y4-'^{9,<j))] function. We have 



{9, (j))=kl + k'^- Qklkl = sin" 9 cos 40 . (2) 

In the gap function (1), the amplitudes of s- and g- 
components are assumed to be equal. Thus 4 (and only 
4) point nodes at 6* = 7r/2 and = 0, 7r/2, tt, 37r/2 are 
realized, see the middle panel of Fig. I. This is exactly 
what has been observed experimentally.® 
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FIG. 1. The xy-plane {9 — n/2) polar plots of the s-|-g 
wave gap |Ak| for various tuning parameter x. 



On the other hand, there is no symmetry reason for the 
constraint of equal amplitudes of s and g. More generally, 
the s-|-g gap function is: 

Ak = A[(l -x) - a; sin"* 61 cos 4(/)] = Af{9,(j)) (3) 

with a tuning parameter x characterizing the weight of 
g-wave component. Obviously Eq. (1) corresponds to the 
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special case x=l/2. If a; < 1/2, s-wave is dominant and 
the nodes will bo removed; while if a; > 1/2, g-wave has 
strong weight and will contribute 8 line nodes. The three 
different cases have been shown in a^y-plane in Fig. I. 

A natural question then arises how to understand the 
origin of the above s+g hybrid pairing. So far, a mi- 
croscopic theory for the pairing potential is not available 
which might be complicated due to the complex Fermi 
surface of borocarbides^^ and the possibility of strongly 
anisotropic Coulomb interactions. As a first step, how- 
ever, it is necessary to investigate phenomenologically 
how the s+g state can be realized by constructing an 
appropriate pairing potential. This is the topic of the 
present work. Within BCS weak-coupling theory we will 
show how a stable coexistence of s- and g-wave and in 
particular how the fine tuning a: = 1/2 (s-l-g) can be re- 
alized almost independent of temperature below Tc- 

A similar issue has been addressed by Lee and Choi^^ 
in their theory to explain Raman scattering data, but no 
complete analysis was presented. First, they do not con- 
sider the ^-dependence of the g-wave part, i.e., implicitly 
assume a gap function with cylindrical symmetry. Sec- 
ond, they use gap models with strong g-component. As 
shown before in Fig. 1, these will exhibit 8 line nodes, 
which are inconsistent with the experimental results^ 
which have revealed 4 point nodes. 



II. PAIRING INTERACTION AND BCS THEORY 
FOR THE S+G WAVE STATE 

In view of the orthogonality of s- and g-wave functions 
in Eq. (3), one may naturally express the pairing poten- 
tial as sum of two separable parts whose weight is given 
by two parameters Vg and Vg-}^ 



(4) 



Here the unprimed and primed angles correspond to k 
and k' respectively. Each term in Eq. (4) is separable 
with respect to wavevectors k and k'. 

Here we try to propose a pairing potential, similar to 
that used in Ref. 17 for the cylindrical gap, by adding 
another mixing term Vsg: 



Vkk' 



-[Vs + Vgi^^'\9,^)i:^^H0',^') 

+F,,(VW(^,<^) + V'''^(^', </>'))] • 



(5) 



The reason for the choice of this type of interaction as 
well as its derivation will be discussed later. We simply 
adopt it as the model pairing potential for the moment. 
The standard BCS gap equation is then written as 



A, 



2£k' 



tanh(/3£;k'/2) , 



(6) 



where /3 = l/iksT) and E)^ = v'(ek - A*)^ + is the 
quasiparticle spectrum, is the free electron dispersion. 



and /i is the chemical potential. It is easy to check that 
the gap function (3) is a self-consistent solution of Eq. (6) 
with inclusion of the pairing interaction (5), if the gap 
amplitude A and tuning parameter x satisfy the following 
self-consistent equations: 

l-a. = KE%^t^^MW2) 

+Vsg V V^'^ (^k, 0k)^^%^ tanh(/3£k/2) , (7) 
k 2£k 

-x = VgJ2 ^^'^ (^k, tanh(/3i?k/2) 



2Ek 



2Eu 



tanh{l3Ek/2) 



(8) 



where the angles have been indexed by their correspond- 
ing wavevector for clarity. Replacing the summation by 
integration according to 



V 47r J_n^^ J 



where Uujd is an energy cut-off to enforce the constraint 
kk — /'I < huj£) (ud- Dcbye frequency for phonon- 
mcdiated SC), and A''(0) is the density of states at zero 
energy for the spectrum £k — we may obtain the fol- 
lowing equations: 



l-X=^iVsIl+V,gl2) , 
-X=^{Vgl2 + Vsgh). 



(9) 
(10) 



Above, V, = N{0)V„ Vg = N{0)Vg, V^g = N{0)Vsg are 
redefined dimensionless interaction constants, and the in- 
tegrals Ii^2 are written as follows 



^^/tanh(^^^7V2) 

^£2 + A2/2 



Jo J + A2/2 ^ ^ 

where we use the abbreviated symbols / and ip^^\ and 
Uwd has been taken as the energy unit. 



III. NUMERICAL RESULTS 

We first consider Vgg = 0, i.e., we assume the pairing 
interaction (4). It was found that one or two s-|-g solu- 
tions (i.e., A > 0, < a; < 1) may appear when is a 
few times greater than Vg . On the other hand, it is easy 
to check that the pure s-wave (a; = 0) and pure g-wave 
(a; = 1) are always trivial solutions. Thus one needs to 
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compare their free energies to find the stable solution. In 
unit of N{0){huJD)^ the free energy is given by 

F= — C de I + + I ln(l + g-'^V^'+^'-f^ 

Jo J P 

Detailed calculation shows that the s+g mixed state is 
unstable in most of the parameter space. As an exam- 
ple, we have shown in Fig. 2 all the solutions and their 
relative energies as functions of Vg for fixed Vg = 0.2 and 
T = 0. Two s+g solutions may be present, as shown by 
the dotted and dashed lines. But compared to the pure 
s- and/or g-wave solutions, they are found to be energet- 
ically unfavorable, sec the lowest panel in Fig. 2}^ Thus 
the s+g mixture seems very unlikely under the pairing 
interaction which is the sum of two separable parts (4). 
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FIG. 2. The values A, x, and relative energies AE as func- 
tions of Vg for Vs = 0.2 and T = 0. The solid and dot-dashed 
lines are for pure g- and s-wave solutions, respectively. The 
dotted and dashed lines are for two possible s+g mixed solu- 
tions. The energy of the pure s-wave solution is taken as the 
reference point in the lowest panel. Energy unit for A is Hlod 
and for AE is iV(0)(awD)^ 

Once Vsg > 0, the above situation changes substan- 
tially. In this case the pure s- and g-wave are no longer 



solutions of Eqs. (9) and (10). Furthermore only one sta- 
ble s+g solution is present. We explain the details in the 
following. 

First we consider a special case: Vg = Vg = —Vsg. 
Then, by adding Eqs. (9) and (10) one can immediately 
obtain the solution with x = 1/2 independent of tem- 
perature. This means that irrelevant of temperature s- 
and g-wave coexists always with equal amplitudes. This 
result is obvious because in this case the pairing inter- 
action (5) can be simply factorized again into the form 
f{e, (p) * f{e', (j)') with fixed a; = 1/2. Then only the gap 
amplitude A(T) is left. It decreases gradually with T 
and vanishes at the transition temperature T^, as shown 
by the dashed line in the upper panel of Fig. 3 where 
Vs = 0.2 is used. 
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FIG. 3. The order parameter A and tuning parameter x as 
functions of temperature T. A and T are in units of Hljo- 
The solid lines are for V = 0.2, Vg = 0.1, Vg = -0.22, and 
the dashed lines are for Vs = Vg = —Vsg = 0.2 corresponding 
to the special case of T- independent a; = 1/2 s+g wave. 

Generally, the three interaction parameters have differ- 
ent absolute values. Various situations can be described 
by tuning these model parameters. We use the follow- 
ing strategy to obtain reasonable values: First, we fix 
the value of Vs, e.g. 0.2 throughout the work which sets 
the overall scale for Tc- Then we assume a value for Vg 
smaller than Vs and tune Vsg to realize the s+g state with 
the constraint: a; = l/2atr = 0. Experimentally, the 
detection of nodes by the field angular-dependent ther- 
mal conductivity is applicable only at very low tempera- 
tures, i.e., T <C Tc. Thus the measurement actually pro- 
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vidcs this constraint on the gap function. In this way, the 
value Vgg may be determined for each given Vg. We have 
obtained a nearly linear relation between Vgg and Vgi 
Vsg ~ —0.24 + 0.2Vg. With interaction parameters fixed 
we can now study the intriguing issue how the x = 1/2 
fine tuning s+g state at T = evolves with temperature. 
One would expect that not only the gap amplitude A, 
but also the tuning parameter x will change with tem- 
perature. If it would decrease, the node points would 
cease to exist and a gap would open with increasing tem- 
perature. In principle this is indeed observed. As an 
example, we take Vg ^ 0.1 < Vg. Then Vsg ~ —0.22 is 
obtained to realize the a; = 1/2 s+g solution at T = 0. 
Under these interaction parameters, A(T) and x{T) are 
calculated self-consistently from Eqs. (9) and (10). The 
results are shown by the solid lines in Fig. 3. It is in- 
teresting to see that x varies with T monotonically and 
very slowly. In the current example with Vg <Vs, x be- 
comes less than 1/2 at finite T. However, the deviation 
from 1/2 is less than 1% even at T = Tc. This means 
that very strong anisotropic gap with Amin/'^max <10~^ 
is present in the whole superconducting region describing 
essentially still a gap with point nodes. Thus we conclude 
that s+g pairing with nodal excitations is a robust solu- 
tion for all temperatures and should not be considered as 
accidental. We also mention that x becomes larger than 
1/2 at finite T for > Vg (not shown), but again only 
a small deviation is obtained. 



IV. DISCUSSION AND CONCLUSION 

In the above section, we have obtained the stable s+g- 
wave based on the phenomenological pairing interaction 
(5). Also we have realized the a; = 1/2 hybrid state which 
has point nodes by appropriate choice of the interaction 
parameters, and most importantly have proven its ro- 
bustness with temperature variation. A microscopic jus- 
tification for our phenomenological scenario is yet to be 
investigated. On the other hand, we should return back 
to understand the form of the pairing potential Eq. (5) 
in more detail. This is discussed in the following. 

For a system with full translation symmetry, the pair- 
ing potential T^ck' can only be a function of k — k'. Fur- 
ther it may be considered to be a function of the angle 
between the two wavevectors in view of |k|, |k'| ~ kp (the 
Fermi wave vector). Then Vkk' = V(k ■ k'), with k • k' 
denoting the cosine of the angle between k and k', can 
be expanded in terms of Legendre polynomials. And by 
use of the spherical harmonic addition theorem, it can be 
finally written into 

oo I 

V{k-k')=Y.^i J2 Yim{9A)Yil{e',^') (13) 

1=0 m=-l 

through spherical harmonics Yimid, 0)- Obviously the Vg 
term corresponds to / = and the Vg term results from 



the sum of (I, m) = (4, 4) and (4, —4). However the cross 
term ^ V^g in Eq. (5) cannot be obtained from Eq. (13). 
This leads us to go beyond the above assumption, i.e., to 
consider Vkk' as a general function which depends on k 
and k' individually. This will be the case when the effect 
of having only discrete lattice translation and rotation 
symmetry is included. Then the pair potential has to 

be expanded in terms of the basis functions V'r {d, (p) of 
the crystal symmetry group (04;^) belonging to a specific 
irreducible representation V of degree I and degeneracy 
index i. The generalized expansion then reads, suppress- 
ing the multiplicity index of F: 

v^^, = Y.^f\k,k')Y,'4^f\e,^)4''\e',ci^r- (14) 

TW i 

For basis functions of different degree /, I' but belonging 
to the same representation F the contribution will gen- 
erally be nonzero. Then for r= Aig and Z = 0, I' = 4 
or vice versa one can naturally obtain the nondiagonal 

contributions V^°^^ = Vgg in symmetry. 

In conclusion, the s+g hybrid gap function has been 
studied within BCS theory. By an appropriate ansatz for 
the pairing interaction and selection of the interaction 
parameters we have shown that the stable coexistence 
of s- and g-wave, as well as the nodal state with equal 
amplitudes of them is possible. In particular, we have 
confirmed that the hybrid state with highly anisotropic 
gap may be robust in the whole superconducting region. 
This provides the semi-phenomenological theory to un- 
derstand the s+g {x = 1/2) gap function proposed by 
Maki et al.^^ Thus the explanation of the experimentally 
observed point nodes in borocarbides Y(Lu)Ni2B2C can 
be achieved in a self-consistent way. Finally we mention 
that similar s+g mixed gap functions are also proposed 
very recently for skutterudite PrOs4Sbi2,^'' and their jus- 
tification should be possible based on a similar analysis 
as presented here. 
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